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ABSTRACT 

Neutrino-neutrino interactions inside core-collapse supernovae may give rise to collective flavor 
oscillations resulting in swap between flavors. These oscillations depend on the initial energy 
spectra, and relative fluxes or relative luminosities of the neutrinos. It has been observed that 
departure from energy equipartition among different flavors can give rise to one or more sharp 
spectral swap over energy, termed as splits. We study the occurrence of splits in the neutrino 
and antineutrino spectra, varying the initial relative fluxes for different models of initial energy 
spectrum, in both normal and inverted hierarchy. These initial relative flux variations give rise to 
several possible split patterns whereas variation over different models of energy spectra give similar 
results. We explore the effect of these spectral splits on the electron fraction, Y e , that governs 
r-process nucleosynthesis inside supernovae. Since spectral splits modify the electron neutrino and 
antineutrino spectra in the region where r-process is postulated to happen, and since the pattern 
of spectral splits depends on the initial conditions of the spectra and the neutrino mass hierarchy, 
we show that the condition Y e < 0.5 required for successful r-process nucleosynthesis will lead to 
constraints on the initial spectral conditions, for a given neutrino mass hierarchy. 
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1 Introduction 



Neutrinos from a core-collapse supernova can play an important role in probing both neutrino 
properties as well as throwing light on the supernova mechanism [1, 2]. Neutrinos emitted during 
the explosion of core-collapse supernovae (SN) pass through very large density variation of matter 
and can undergo MSW resonant flavor conversion [3] which can give useful information on neutrino 
mass hierarchy and the third leptonic mixing angle 6*i 3 . They are also influenced by the shock wave 
formed in the SN core and thus carry information about the explosion mechanism [4, 5, 6, 7, 8, 9]. 

Recently it was realized that a crucial feature in the study of SN neutrinos comes from the 
collective neutrino-neutrino interaction at very high densities of the core and this may change the 
emitted flux of different flavors substantially [10] -[42]. The phenomenology of supernova neutrinos 
with collective effect including Earth matter effect on SN neutrinos [35], prompt SN [31, 32], 
diffuse neutrino background from SN [43], failed SN [44] as well as CP violation in SN neutrinos 
[38] has been studied. For a recent review we refer to [45]. 

It has been shown that effectively the collective evolution of a three-flavor (z/ e , z/ M , u T ) system 
can be treated like a two flavor (v e , v x ) scenerio, where v x can be or v T or a linear combination 
of and v T [33, 40]. The flavor evolution for this system is driven by the effective mass squared 
difference Am 2 and the mixing angle #13. This two flavor scenerio which has been extensively 
studied [20, 25, 26, 29] shows that for inverted hierarchy (IH, Am 2 < 0), above a critical energy 
(split energy E c ), the spectrum in both the electron neutrino (i/ e ) and antineutrino (P e ) sectors 
end up with a complete exchange or swap with v x and V x respectively, this is referred as "spectral 
swap". Recently the studies in [41, 42] analyzed the role of equipartition in energy and variation 
of luminosity and showed the interesting possibility of multiple splits in the supernova neutrino 
spectrum for IH. Single spectral split for Normal Hierarchy (NH, Am 2 > 0) was also reported for 
certain values of luminosities. 

In this paper our goal is to study (a) the impact of spectral splits on the electron fraction 
(Y e ) which is a diagnostic of successful r-process nucleosynthesis in supernova and (b) the inverse 
problem which is to see if it is possible to put any constraint on the luminosities by demanding 
the neutron rich condition on Y e which is required for synthesis of heavy nuclei. 

Though the site for r-process nucleosynthesis is not known definitely, supernovae are considered 
to be excellent candidates for it. One of the criteria for the rapid nucleosynthesis to take place is 
that it has to be in a neutron-rich region. With the two competing beta processes n + u e — > p + e~ 
and p + v e — > n + e + occurring in the hot bubble and neutrino driven wind region, the minimal 
condition is that the electron fraction, Y e < 0.5. A more realistic constraint may be Y e < 0.45. 

Effect of MSW neutrino oscillations on r-process nucleosynthesis were considered in [46], and 
subsequently in a large number of papers. Since r-process is expected to take place in the neutrino 
driven wind deep inside the supernova, the matter density in these regions obviously are high. 
Therefore, in order to have MSW effect one requires mass squared difference Am 2 ~ 1 eV 2 or 
higher. Such large values of Am 2 are possible only if one allows for sterile neutrinos. For only 
active neutrinos, the mass squared differences are of the order of 10 -3 eV 2 (atmospheric) and 10 -5 
eV 2 (solar) only, and the MSW resonance for these are reached at distances far beyond the r-process 
region. Therefore with only active neutrinos, one expects no effect of MSW oscillations on r-process 
nucleosynthesis. However, one will have collective effects driven by these active neutrino mass 
squared differences and these collective flavor oscillations happen very close to the neutrinosphere 
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- close enough to impact r-process nucleosynthesis. The effect of the collective oscillations on the 
possibility of getting n-rich region in the hot bubble was studied in [17]. The problem of r-process 
in the neutrino- wind region with different evolution scenarios was looked into in [47]. With an 
improved understanding of the spectral splits in the collective oscillations [41, 42], we re-examine 
the problem of SN r-process using the two flavor basis as mentioned above. In particular, recent 
papers have shown that (multiple) spectral splits could happen in either or both the neutrino and 
antineutrino channels, depending on the initial spectral conditions. Since spectral splits can change 
the neutrino and antineutrino spectra and hence the value of Y e , and since these spectral split 
patterns depend on the initial conditions of the unoscillated spectra, r-process nucleosynthesis can 
in principle be used to put constraints on the initial spectral conditions. We perform this exercise as 
a proof of principle, in a simplified framework, in order to illustrate that such constraints indeed 
exist. We consider different models for the initial neutrino spectrum with different hierarchies 
among average energies and departures from energy equipartition among flavors. We demonstrate 
the occurrence of spectral splits due to these variations on the probabilities and the fluxes for both 
hierarchies. We use the constraint that the environment should be neutrino rich for successful r- 
process and delineate the resulting constraints that are obtained on the fractional luminosities. 

The paper is organized as follows. In Section 2, we outline the neutrino-neutrino interaction in 
the context of the supernova problem. First the two flavor evolution equations are presented. Next 
we vary the fractional luminosities for different initial energy spectrum model and study the the 
effect of collective oscillations, especially the effect of the spectral splits on the emitted neutrino 
flux. In Section 3, we discuss how the evolution of the electron fraction Y e is affected by collective 
oscillations. The value of Y e determines the possibility of having r-process nucleosynthesis in the 
supernova environment. We impose the constraint that the environment should be neutron rich 
and study the initial relative fluxes or relative luminosities which are allowed under this constraint. 
Finally Section 4 makes some concluding remarks. 

2 Two flavor neutrino-neutrino interaction 
and Supernova 

Close to the neutrinosphere, due to the large neutrino density, neutrinos form a background to 
themselves. This neutrino-neutrino interaction effect is nonlinear and can give rise to collective 
flavor transition of neutrinos and antineutrinos. In this section we discuss the flavor evolution 
equations, different models of supernova neutrino spectrum and the impact of the collective effects 
on the probabilities and the fluxes for both NH and IH. 

2.1 Two flavor evolution equations of SN Neutrinos 

Due to the large neutrino density inside the neutrinosphere, the neutrino-neutrino interactions 
lead to coherent oscillations of neutrinos of different energies with some average frequency, giving 
rise to synchronized oscillations. However, there is no effective flavor conversion due to these 
synchronized oscillations as the effective mixing angle is highly suppressed due to the large MSW 
potential in the region close to the neutrinosphere. With the neutrino density decreasing out- 
ward, bipolar oscillations begin to take place. These oscillations can lead to complete or partial 
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swapping (spectral split) of the v e (u e ) and v x (u x ) spectra depending on their initial luminosities 
and average energies. Finally, after a few hundred kilometers, the neutrino-neutrino interactions 
become negligible and it is the MSW transitions which dominate. 

As shown in [33, 40] the collective effect due to neutrino-neutrino interaction effectively involves 
only two flavors of neutrinos (u e , u x ), while the other flavor (u y ) does not evolve under this collective 
potential. Here v x is a linear combination of v T and v y is the orthogonal combination to v x . 
The only way v y can effect the final neutrino spectrum is by MSW transition which happens at a 
larger radius of about 10 4-5 km, well beyond the collective region which is within a few hundred 
km from the center of the exploding star. Therefore, the effective evolution of the neutrinos is 
very well described by the two flavor formalism. 1 

The evolution equations in the two-family Bloch vector notation for the polarization vectors 
of the neutrino (P) and antineutrino (P ) sector are, 

P = Px(uB-Az-/iD), (1) 



P = P' x (- W B - Az - /iD) , (2) 

where the terms involving u, A and ji are the ones having the vacuum, matter and neutrino- 
neutrino interaction effects and the frequencies are represented by 

B = (— sin 20, 0, cos 29) T , u = — — , (3) 

2E 

z = (0,0,l) T , \ = V2G F N e , (4) 
J dE{nP - nP') , n = V2G F (N Ue + N Ux + N Pe + N Px ) , (5) 



D 



{N Ve + N Ux +N De +N Dx 



respectively. 

D can be defined in terms of the global polarization vectors J and J i,e D = J — J . 
Where 



JdEnP , .1 v . v ' . v , /r//w/P. 



(N Ve +N Vx +Nv e +N& X ) J ' (N Ve +N Vx +N Be +Np x ) 

As usual, 9 and Am 2 are the mixing angle and mass squared difference respectively. In what 
follows 9 e ff is taken as 1CT 5 and |Am 2 | = (Am^l = \m\ — m\\ =3 x eV 2 . N^s represent the 
total effective number density of the ath species. 

N a = JdEn a , (6) 



1 The MSW [i-T resonance can occur inside the supernova under certain conditions [54]. Unlike the standard 
MSW resonances involving the electron flavor, the \i-t resonance can occur in the inner supernova layers. The 
effect of jU-r neutrino refraction and collective three-flavor transformations in supernovae has been considered in 
[55]. We have ignored these effects in this work. 
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where, 



n = n Ue +n Ux , n = n Pe + n Px , (7) 

n Q 's are the effective number density per unit energy for the a'th species of neutrino and can be 
expressed as [20] 

where L a and (E a ) are the luminosity and average energy for the ath (anti)neutrino species, R a 
is the neutrinosphere radius. The initial flux of the ath species at the neutrinosphere is given by 
whereas the initial energy distribution is represented by ty a (E). D(r) denotes the geometrical 
function in the "single angle approximation" [20] can be expressed as, 2 




D(r) = - U-\ l- — • (9) 



The matter effect is removed from the evolution equations as the equations are considered in a 
frame rotating with angular velocity -Az [19]. 3 In such a frame all the physical observable remain 
the same. Thus the evolution equations are 

P = Px(wB-/iD), (10) 



P = P' x (-wB - fiD) . (11) 

These are nonlinear coupled equations (due to the 2nd term containing D) and have to be solved 
numerically. It is evident from the evolution Eqs. (10) and (11) that there are two relevant frequen- 
cies, the usual vacuum frequency (uu) and the neutrino-neutrino interaction strength parameter 

2 It is well known that the current-current nature of the weak interaction introduces an angular factor in the 
Hamiltonian because of which the neutrinos from the SN core traveling along different trajectories encounter 
different neutrino-neutrino interaction potential. These 'multi angle' effects may give rise to kincmatical decoherence 
[22] which in turn can wash out the collective features described above. But for spherically symmetric cases 
"single angle" approximation i,e neutrino-neutrino interactions averaged along a single trajectory seems to be a 
fine approximation as the 'multi angle' decoherence in this case is rather weak against the collective features [23]. 
We are aware that the extent of kinetic decoherence is initial flux condition dependent and hence some of the initial 
fluxes that we consider might lead to kinetic decoherence. In [42] the different possible split patterns with varying 
neutrino luminosities were studied using the "single angle approximation" . In this work we study the multiple split 
patterns for different models of initial supernova neutrino spectrum in the same spirit as [42] . Multi angle study of 
these possible spectral patterns for different initial supernova neutrino spectrum might give interesting results like 
the effect of decoherence, but this is beyond the scope of our work. 

3 Note that though the matter potential is mostly rotated away, it may affect the evolution by delaying the 
collective effect [21] or by some early decoherence [36] or even modifying very low energy (order of 0.1 MeV) 
split features [29, 30, 2]. These early effects have very little impact on the over all split patterns at the end of 
the collective region (400 Km) and the low energy split features below 1 MeV are negligible compared to the 
total spectra. Moreover the matter potential can be accounted for by choosing a matter suppressed hence small 
effective mixing angle, as we have chosen a e ff = I0~ 5 [23]. So in the subsequent discussions we neglect the above 
mentioned roles of the matter term and work with a very small e ff to compensate the matter term. We have 
explicitly checked that the inclusion of the matter term does not change our results. 
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([/,). For our chosen Am 2 the vacuum frequency is 



Am 2 30 , 

^ = ^ - = 7TTTTTT T r , km . (12) 

2£ (AE/MeV) v ; 

The usual SN neutrino energy considered is in the range of to 50 MeV, as the SN neutrino flux 
beyond 50 MeV is very small. Hence we use neutrino energy upto 50 MeV for our calculation. 
The other frequency (/i) representing neutrino-neutrino interaction is nontrivial, and is given 

by 



/i 



= V2G F (N Ve + N Ux + N Pe + N Px ) . (13) 



The Eqs. (6) to (9) imply that contribution of the a-th species to /i is dependent on radial distance 
(r), neutrinosphere radius (R a ), initial flux (y^fy) and initial energy distribution (^ a (E)). 
In our analysis, neutrinosphere radius is taken as 10 km whereas other inputs like initial flux 
and energy distribution depend on the choice of initial neutrino spectrum model. We analyze the 
evolution for several neutrino spectrum model. 



2.2 Models of Initial Neutrino Spectrum 

In a core-collapse SN, the gravitational binding energy (about a few times 10 53 erg) is converted 
to neutrinos and antineutrinos with energies of the order of 10 MeV and gets emitted in the 
subsequent ~ 10 sec. Initially a neutronization burst comes out consisting of pure z/ e s but with only 
a very small fraction of the total energy and after that the thermal neutrinos and antineutrinos of 
all three flavors are emitted. For the thermal neutrinos the initial energy distribution is expected to 
be Fermi-Dirac (FD), but the results of several simulations [48, 49, 50] found that the distribution 
must be close to pinched thermal spectra [50] i.e. with a deficit on the high energy side compared 
to FD. Fermi-Dirac(FD) distribution in energy implies 

- -^sff ' (14) 

and for a choice of average energies of different flavors 

(E Ve ) — 10 MeV, (E Pe ) = 15 MeV, (E Vx ) — (E Px ) — 24 MeV, (15) 
the inverse temperature parameters are [29] 

p Ve = 0.315 MeV" 1 , (3 Pe = 0.210 MeV^ 1 , (3 Ux = (3 Px = OTSlMeV" 1 . (16) 
Whereas the pinched spectra for different simulations are parameterized as [50] 

(E a ) is the average energy of u a , and ( a is the pinching parameter. 
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The effective number density for the ath species per unit energy is given by 



"«c< £ > = S|<t>w < 18 > 

For a specific choice of ^ a (E) the initial flux (<fi a = -^y) for the crth flavor need to be specified 
and are very crucial input parameters in our study. Supernova models tell us that almost all the 
gravitational energy released in core collapse supernovae comes out as uus of all flavors. Only 
one or two percent of it goes into the explosion and the electromagnetic radiation emitted in all 
wavelengths. The total luminosity scales as L(t) = L (e _ '/ r /r) but for a first study we take a 
time-averaged value for it as done in [29] and [42]. One can of course look at the problem for 
specific instants of time by changing the total luminosity, early times having larger values. 

The total SN binding energy released (Eb = 3 x 10 53 erg) is related to the individual flavor 
luminosities by 

L Ve +L Pe +4L Ux = ^ , (19) 

T 

assuming no distinction between v x and v x . We also assume a time- independent constant lumi- 
nosity over the time r. We take r = 10 seconds. Thus the initial fluxes of different flavors get 
constrained by 

<Pl(E Ue ) + <Pl(E Pe ) + m x (E Vx ) = 3 x 10 52 . (20) 
If we denote the ratio between the initial fluxes of different flavors by 

< : A ■ = C : ■■ 1 ' (21) 
where (f> r Ue , (f> r Pe are positive numbers, then Eq. (20) can be written as 

tiM e (E„ e ) + ^(E,.) + 4(E Ux }) = 3 x 10 52 . (22) 

0° 

Note that (j) r v = , <p r 9 = -nr are basically initial relative fluxes. Thus different choices of <p r v 

and 0p e imply different relative luminosities or relative fluxes. 

Four representative sets for the energy spectra (in terms of (E a ) and the pinching factor £„) and 
flux ratios usually discussed in literature, are given in Table 1. One simulation by the Lawrence 
Livermore group (LL) [48] and two different simulations by the Garching group (Gl, G2) [50] are 
presented. Recently [41] used another set of "plausible" flux parameters giving rise to multiple 
splits in the neutrino spectra is also given. We call this 'G3'. For the LL spectra we use the FD 
distribution for \1/ given in Eq. (14). The ft a for LL are given in Eq. (16). For Gl, G2 and G3 
spectra we use the pinched spectrum defined in Eq. (17). We assume ( Ux — („ x = 4 and ( Ue =(p e = 3 
for Gl and G2. For G3 all ( a = 3. 

Note that the LL simulation obtained a large hierarchy (E„ e ) < (E Pe ) < (E Vx ) pa (E Px ), and 
an almost complete equipartition of energy among the flavors. The Garching simulations predict 
a smaller hierarchy between the average energies, incomplete equipartition, and increased spectral 
pinching. The differences in the values of these parameters arise from the different physics inputs. 



7 



Model 


(Eu.) 




\Ev x ,v x ) 


€ = 


€ e 






LL 


12 


15 


24 


2.00 




1.60 




Gl 


12 


15 


18 


0.80 




0.80 




G2 


12 


15 


15 


0.50 




0.50 




G3 


12 


15 


18 


0.85 




0.75 





Table 1: The parameters of the used primary neutrino spectra models motivated from SN simula- 
tions of the Garching (Gl, G2) and the Lawrence Livermore (LL) group. We assume ( Ux =(u x = 4 
and ( Ue =(p e = 3 for Gl and G2. For G3 all ( a = 3. For LL we use a pure FD spectrum. 



The equipartition of energy implies 

L Ue = L Ue = L Ux . (23) 

In terms of our notation it means 

So complete equipartition for the Garching simulations would imply flux ratios (Table 2) different 
from the values in Table 1. Recent analyses [42] have shown that the multiple split cases have 
origin in the departure from energy equipartition. 

Actually there is no reason that equipartition should be strictly followed for the energy released 



Model 






<Pk = 


< 


Gl 


1.50 




1.20 




G2 


1.25 




1.00 




G3 


1.50 




1.20 





Table 2: The flux ratios for the Garching models with equipartition of energy 



from a real supernova. In the next subsection we make extensive analysis of this multiple split 
phenomena with varying initial fluxes, which is equivalent to varying <p r Ve and 0p e . 

2.3 Survival Probability and Flux 

As stated above in this subsection we discuss the impact due to the variation of initial relative 
fluxes ((jf v& and (jf De ) on the final spectrum. The final spectrum is calculated at 400 km as collective 
effect is expected to vanish at around 400 km. We also analyze this effect for different models of 
initial neutrino spectrum spectrum LL, Gl and G3. 

In principle the values of <p r Ue and <f) r Pe can lie in a large range. Thus analysing this variation 
would require study in a wide range of the </£ e -0p e parameter space. Instead we consider the 
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suggestion [51] that the uncertainty in the relative luminosities of different flavors must be in the 
range 

-<^l<2 ; -<^<2. (25) 
2 L Vx 2 L Vx 

These limits in turn will put a constraint on the parameters 0^ and 0p e 

2 jej £ < £ 2 m • 2 je-) 2 ** 5 2 m ■ (26) 

In Table 3 we present the lower limits (11) and upper limits (ul) of the initial relative fluxes for 
different spectrum models LL, Gl and G3. 



Model 


(E„ e ) 










0c 




LL 


10 


15 


24 


1.20 


4.80 


0.80 


3.2 


Gl 


12 


15 


18 


0.75 


3.00 


0.60 


2.4 


G3 


12 


15 


18 


0.75 


3.00 


0.60 


2.4 



Table 3: The average energies, upper limits (ul) and lower limits (11) of the initial relative flux for 
the models used. 



To compare different flux models and study more of the parameter space we vary <p r Ue and <jf Ds 
in the range [0.5,5.0] and [0.5,3.5] respectively, for all the models. We find that varying 4> r Ve and 
0p e give rise to different possibilities of final spectra as discussed in [42]. In addition to that we 
check it for different initial spectrum models. Here it is notable that usually the initial spectrum 
models come with a fixed value of <p r Ue and (jf 9e (see Table 1) but the main idea in this analysis is 
about varying <p r Ue and 0p e . So here by initial spectrum models (like LL, Gl, G3) we mean the 
energy dependence (( Ua ) and neutrino average energies ((E Ua )) of the models. In what follows, we 
will see that for the Inverted Hierarchy (IH) the final spectrum is very sensitive to the values of 
c/£ e , 0p e and the model of initial spectrum. Whereas for Normal Hierarchy (NH), the results are 
less dependent on these quantities. We will discuss the reasons for this behavior. 

2.3.1 Probability and Flux: NH 

As already discussed in [41, 42], large flux (luminosity) of v x can induce simultaneous swap in 
both neutrino and antineutrino sector for NH. In these cases initially the system is in an unstable 
equilibrium. As it evolves, it partially swaps the flavor in both neutrino as well as antineutrinos, 
to end up in a stable state. We further study this over different spectrum models and initial 
relative fluxes. For each of the different models we vary <f) r Ue and <p r Pe in the range [0.5,5.0] and 
[0.5,3.5] respectively. We find that for several choices of {4>^ e ,4> r p e ) there is simultaneous swap in 
both neutrino and antineutrino spectrum and this swap may generate prominent split in the final 
spectrum. For a specific model these split energies (E c ) may vary from low to high energies, 
depending on the value of (</£ e ,</>p e )- 
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Figure 1: Survival probability for G3 spectrum in NH with different {4> r v ,cjf 9 ). 
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Figure 2: Survival probability for LL and G3 spectrum in NH with different 



Independent of the choice of spectrum models, these split features are seen for low values of 
6p ), which implies large flux of u x compared to other flavors [42]. As the values of {4>^ e ,4> r 9e ) 
increase, the split energy (E c ) also increases, and close to the equipartition point the split energy 
tends to infinity. 

As an example see Figure 1 where survival probabilities are plotted for the spectrum G3. The 
left panel is for neutrino and right one for antineutrino. For a low value of (f) r Ve (0.5), (ff p< _ is increased 
from 0.6 to 2.4 for both neutrino and antineutrino. ^From the figures it is evident that with this 
increment, split energies (E c ) also increase. For the same combination of <p r Ue and 0£ e the split 
energy is higher for neutrinos than for antineutrinos. We find these features are same for other 
spectrum models too. 

In Figure 2 we have shown the above mentioned features for LL and G3. The left panels are for 
LL and the right ones are for G3. The red straight lines are for neutrino and the blue dashed lines 
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Figure 3: Flux for the different neutrino species for G3 with NH in arbitrary units (a.u.); WOC stands for 
"Without Collective" effects and WC for "With Collective" effects. 

are for antineutrino. In the top panels the initial relative spectrum (0£ ) are chosen to be close 
to the equipartition point for both the models and for these values there is no split whereas for the 
lower panels ,4> r p ) are smaller and these panels show split for both neutrino and antineutrino 
sectors. 

The flux corresponding to G3 and (0.5,0.6) are plotted in Figure 3. Left and right panels in this 
figure are respectively for neutrinos and antineutrinos. The upper panels are without collective 
effects (initial flux) and the lower ones are with collective effect (flux beyond collective region). 
The black lines in the panels are for electron type whereas the red dashed lines are for x-type. 
Clearly the lower panels show swap in both neutrino and antineutrino sector. The swap in both 
sectors are partial, that is, a part of the spectra below the "split energy" remains same. The 
antineutrino split feature is not clearly visible since E c for them is low and the v e and u x fluxes are 
very close to each other at these energies. The probability plots for (0.5,0.6) in Figure 1 also show 
the low split energy for antineutrino. Note that the swap for neutrino spectra happens at a higher 
energy compared to the antineutrino spectra, this feature is also consistent with the probability 
plots in Figure 1. 

2.3.2 Probability and Flux: IH 

Probability and flux in the IH is much more complex and interesting than NH. Here also we 
vary the initial relative flux for different spectrum models and find wide variation of the final 
spectrum depending on the choice of These variations in spectrum with initial relative 

flux ,<fip) have been attributed to meeting the instability condition of the initial system and 
adiabaticity violation [41] as well as to change of global initial condition with luminosity variation 
and minimization of potential energy [42] . We find that these changes in final spectrum are similar 
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Figure 4: Survival Probability for LL and G3 spectrum in IH with different , 4>p )• 



for different choices of spectrum models. 

As discussed in [42] , the different spectral features arise from the initial conditions, which may 
or may not lead the system to swap to minimize "potential energy". We also find that in some 
cases the multiple swaps actually do take place but the swaps are so close that they can not be 
resolved numerically [41] and thus appear as if the swap or split features are absent. We find five 
spectral split patterns as mentioned in [42]. These five patterns are found for all three models of 
initial energy spectra LL, Gl, G3. These are displayed (for LL and G3) in successive panels from 
top to bottom in Figure 4. 

1. Dual split in both neutrino and antineutrino flux (II, II). 

2. Dual split in neutrino but no split in antineutrino flux (11,0). 

3. One split in both neutrino and antineutrino flux with the split energy of the neutrino higher 
(H) than that of antineutrino (L) split energy (1,1) (H,L). 

4. One split in both neutrino and antineutrino flux with the split energy of the neutrino lower 
(L) than that of antineutrino (H) split energy (1,1) (L,H). 

5. No split in neutrino but dual split in antineutrino flux (0,11). 

Apart from these five patterns we find a sixth possible pattern in which neither neutrino nor 
antineutrino show any swap in the spectrum. We call this (0,0). For this pattern (0,0) the effect 
of neutrino-neutrino interaction on both neutrino and anti neutrino flux is undetectable. 
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Figure 5: Survival Probability for LL (3.4, 0.5) and G3 (2.3, 0.8) spectrum in IH . 



The physical reasoning behind the patterns in the top five panels are well explained [42] from 
the idea of potential energy minimization. Our analysis shows that in some sense all the different 
spectrum models are in the same footing as all of them give rise to similar split patterns with 
the change of initial relative flux or relative luminosity. As explained in [41] the basic feature is 
that there are multiple swaps or splits in both the neutrino and antineutrino sector but the swaps 
may disappear depending on the adiabaticity violation or it may be numerically unresolvable. split 
patterns Consider the new pattern, described in the lowest panel of Figure 4, where it seems that 
there is no swap in both neutrino and antineutrino sectors. When we study these cases carefully 
(Figure 5, left panel LL (3.4,0.5) and right panel G3 (2.3,0.8)) we find that they also show changes 
in survival probability similar to the other patterns. But the swaps here are incomplete and 
numerically undetectable. While in Figure 4 the change in probability for this case is visually 
unresolvable for all practical purposes, in Figure 5 it is visible, as we have increased the resolution. 

For the fluxes we just give one example of the case (II, II) in Figure 6. Here we plotted the G3 
neutrino spectrum in the left panel and the G3 antineutrino spectrum in right one. In both panels, 
the solid sky blue lines are for electron type without collective effect (WOC) and the solid red lines 
are for v x without collective effect (WOC). For the spectrum with collective effects (WC) dashed 
black lines are for electron type whereas dot-dashed blue lines are for v x . Here the spectrum model 
used is G3 and the initial relative fluxes are (1.1,0.8). We can see prominent dual split pattern in 
this flux figure as expected from the upper right panel of Figure 4. 

Thus, depending upon the choice of initial relative fluxes (</>y e ,</>p e ) the spectra can have different 
patterns, especially for IH. The possible values of (0^ e ,0p e ) can be in a wide range. Even if 
one assumes a factor-of-two-uncertainty in the relative luminosity [51], there can be considerable 
variations in the final flux characteristics. 

We study the variation in spectral split features over the 0^ e -0p e plane for LL, Gl, G3 and 
found a pattern showing different kind of spectral splits at different (f> r Uf ,-4>p e region. In Figure 
7 we show this in the (ff Ue -(f) r Ue plane. The plane is divided into zones by the values of the global 
polarization vectors J z , J z and D z . The black dashed line divides the plane into zones with D z > 
and D z < 0. The purple long dashed corresponds to J z = and demarcates the area which has J z 
positive and negative. The blue thick dashed is for J z = 0. These lines therefore divide the </£ e -</>£ e 



13 




Energy (MeV) Energy (MeV) 



Figure 6: The neutrino and antineutrino fluxes in arbitrary units (a.u.) for the model G3 with the relative 
luminosities (1.1,0.8), with and without collective effect for IH. 

plane into 6 zones. The split patterns observed in the different zones are shown on the plane. 
The global polarization vectors J z , J z and D z , define the "phase transitions" across different split 
patterns. It should be noted that the global polarization vectors were initially in the z direction 
hence the sign changes of their z components mark the stability of the system and spectral splits 
[42]. Figure 7 shows that 

1. (1,1) (H,L) patterns are for J z > 0, J z > and D z > 0, 

2. (1,1) (L,H) for J z > 0, J z > and D z < 0, 

3. (a) (II, II) patterns are seen in the J z < 0, J z < region, and also in the J z > 0, J z < 

region, 

(b) (0,0) appear mostly in J z < with a very few occurrence in J z < 0, 

(c) (11,0) is the most dominant pattern in the J z < region, although it can appear in the 
J z > 0, J z < region, 

(d) (0,11) pattern occurs only in J z < 0. 

Thus for the double split patterns described in point 3 above, the so called "phase transition" lines 
seems inconclusive. [42] established that the J z < or J z < region (i.e, the rectangle covering 
the zone with J z < alongwith the rectangle covering the zone with J z < ) lead to (II, II) 
pattern when the adiabaticity is increased artificially. Hence due to the incomplete adiabaticity in 
the actual case some of the (II, II) patterns appear as (11,0), (0,11), (0,0) pattern. So with complete 
adiabaticity one will not see any of the (11,0), (0,11), (0,0) patterns in Figure 7. For example see 
the discussion regarding the (0,0) pattern in context of Figure 4. 

With so many possible patterns it will be really difficult to predict the initial relative neutrino 
fluxes, the energy distribution model and the extent of collective neutrino effect, even for a future 
galactic supernova event. In the next section we discuss to what extent one can constrain the 
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have different split patterns. The abbreviations used are explained in section 2.3.2. 



luminosities by demanding a neutron rich condition by the end of the collective region, required 
for successful r-process. 

3 Neutrino Fluxes and r-Process Nucleosynthesis 

In this section we discuss the effect of the flux of neutrinos radiated out in core collapse super- 
novae on the electron fraction and discuss the possibility of getting allowed regions for r-process 
nucleosynthesis and the resulting constraints on relative luminosities. As most simulations of 
core-collapse supernovae do not lead to explosions, there are uncertainties in the understanding 
of the late stage of the SN shock propagation. But the generally accepted scenario supported 
by simulations is that for core collapse supernovae starting with iron cores the shock wave gets 
initially stalled due to loss of energy through nuclear dissociation and then over timescale of a 
second, gets revived by the energy deposited by neutrinos radiating out, the so-called late-time 
neutrino heating mechanism leading to the delayed core collapse supernova [52]. This results in 
the development of a low density "hot bubble" region just behind the SN shock. Normally the hot 
bubble regions are taken between the infalling neutron star radius and the forward shock, that 
is, up to 30-40 km initially. The huge flux of neutrinos emitted from the proto-neutron star leads 
to the "neutrino-driven wind" which remains active for about 10 seconds after the core bounce. 
This creates neutron-rich regions of high entropy which are conducive to the development of the r- 
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Figure 8: The electron fraction Y e for both hierarchies as a function of 'r', the distance of the region from the 
center of the core. The spectrum model used is G3. 

process. Different delayed core-collapse SN calculations give rise to different values for the entropy 
per baryon leading to confficting conclusions about the r-process. However the //-driven wind is 
still considered to be one of the most probable sites for the r-process [52]. The //-driven wind 
models consider the r-process site to be at a few hundred kilometers (within 1000 km) [52]. Since, 
as discussed before, this is also the region where collective oscillations are active, it is expected 
that r-process will get affected. In what follows, we study this effect in the //-driven wind region 
and numerically check if criteria of successful r-process can be used to constrain initial neutrino 
flux parameters. 

The criteria for r-process on which we focus here is the the electron fraction Y e , i.e., the 
number of electrons (equal to the number of protons, due to charge neutrality) per baryon. The 
Y e will depend on the relative strengths of the two reactions - neutrino capture on neutrons and 
antineutrino capture on protons. Therefore, Y e can be expressed as [46] 

Y e = 1/(1 + X PeP /X Uen ) , (27) 

where X Uen and \ PeP are the reaction rates for v e + n — > e~ + p and v e + p — > e + + n respectively. 
Note that these reactions can in principle occur on both free and bound nucleons. However for 
the purpose of this work, we will not consider the reactions on heavy nuclei 4 . Note also that 

4 For a detailed study of the effect of nuclear compositions on Y e we refer to [16]. Our main conclusions come 
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in principle, the inverse reactions also happen inside the supernova and should be considered. 
However, we neglect the inverse reactions here since the matter temperature of the region is small 
compared to the neutrino temperature as one goes away from the neutrinosphere and has very- 
small effect at radius of 30 km and beyond [46, 53]. The reaction rates \ uN (where N = n or p) 
are given as 

^ L v JZ°v vN (E)f v (E)dE 
XvN ~ W JS°Ef v (E)dE ' (28) 

where L v = (fP v (E v ), and f v denotes the neutrino flux. The cross section used are 

°UEv.) - 9.6 x 10- 44 (^) 2 cm 2 , (29) 



- 9.6 x 10- 44 ( Ep Me y nP J™ 2 , (30) 

where A np = 1.293 MeV is the mass difference between neutron and proton. 

Note that the neutrino flux denoted as /„ in Eq. (28) is the flux including collective flavor 
oscillations. The swap between the active neutrinos due to collective effect can change f v and 
hence Y e . To illustrate this better we show the ratio of the reaction rates for v e and v e explicitly 
in terms of the collective oscillation probabilities 

K eP , , C ^ P (E)Pj(r, E)<j>**,, e (E)dE + f °° <r PeP (E)(l - i%(r, E))^ Ux (E)dE 

XuJ r> J^a v . n (E)PsSr,E)^ v ^(E)dE + J^a v . n (E)^ ' 1 ' 

where P£ e (r,E) and P£ e (r,E) are the anti-neutrino and neutrino survival probabilities with col- 
lective oscillations and are calculated numerically as function of radius and energy. The minimal 
condition for the SN environment to become neutron reach is Y e < 0.5 which translate as the 
condition Ap eP /A I/en > 1. 

Let us begin by understanding the impact of collective flavor oscillations on r-process by 
discussing some limiting cases, n the no flavor oscillation limit, i.e. Pp e = P^ e = 1.0 for all 
energies, Eq. (31) reduces to 

(XpeA <P r p J °°(E-A np )^ Pe (E)dE <j>l e ((E - A np f) Pe 

U^Jnoosc <PlJ^{E + A np ym Ue (E)dE < e ((£ + A np )% 1 > 

Since the average energy of v e is greater than that of v e for all the three SN models that we 
have considered, it is expected that ((E — A np ) 2 ) Pe > ((E + A np ) 2 ) Ue . Therefore under this 
approximation, for 0p e /0^ e > 1, Y e < 0.5 always and r-process can proceed. The condition 
Y e < 0.5 in fact gives Vp./V^ > 0.62 for LL and > 0.88 for G1/G3 . For all values of <p r p J<j) Ue 
greater than this value, r-process can happen while for all values of <^p e /0£ e below this, r-process 
is forbidden. 

from impact of collective oscillations on r-process nucleosynthesis and are not expected to drastically change as a 
result of reactions on bound nucleons. 
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Likewise one could consider the case where we have complete conversion of both neutrinos 
and antineutrinos where Pp e = P^ = for all energies. One can easily show that for this case 
K eP /K e n = ((E - A np ) 2 ) v J((E + A np ) 2 ) Vx . As a result here one always gets Y e > 0.5 as A np is 
positive. Note however that this case never happens in collective oscillations and is therefore not 
realistic. 

Next we consider effect of collective effects on r-process for the realistic case, where the flavor 
conversions are calculated numerically, as outlined in the previous section. In Figure 8 we show 
the electron fraction Y e as a function of the radius (r), for different combinations of (p r Ue and 
(f) r Pe . We have taken the G3 model for the v e and v e spectra. The upper panel is for IH while 
the lower one is for NH. In the upper panel, the green line is for , (fipj of (1.5, 2.0), the 
blue line for (0.8, 0.6), while the maroon line is for (2.3, 0.8). We can note from the figure 
that for IH and (2.3, 0.8) case, there is no flavor conversion due to oscillations for the inverted 
hierarchy. The probability Pp e and P^ for this case was shown in the lowest right-hand panel of 
Figure 4, where we can see that Pp = P£ = 1. Therefore as discussed above, the \y eP /\ Vt , n = 
(0.8/2.3) * ((E - A np ) 2 ) p J((E + A„ p ) 2 )^ e = 0.396, giving Y e ~ 0.72. For the two other cases 
considered with IH we have oscillations due to single and multiple splits. A scan of Figure 7 
reveals that we have double splits in both neutrino and antineutrino channels for the blue line 
with (0.80, 0.60) whereas for the maroon lines of (1.50, 2.00) we have single splits in both neutrino 
and antineutrino channels with the split energy of neutrino lower than that of antineutrino. For 
both these cases we can see very fast oscillations in Y e within the first 200 km, which can be 
attributed to the bipolar collective oscillations. Beyond 300 km the value of Y e approaches a fixed 
value as the neutrino density decreases very fast and the collective effects end. The reason that 
one gets higher values of Y e for both of them after the completion of collective effects compared to 
their values at 30 km can become clear from Eq. (31). For double splits for fixed value of 4>p e /4> Ue 
the contribution from the integrals in denominator in between the split energies is more than the 
corresponding contribution in the numerator making the ratio lower, resulting in higher Y e . On 
the other hand for the single split case of (1.50, 2.00) the low split energies in the denominator for 
neutrinos make the ratio of Eq. (31) lower. 

In the lower panels we assume NH and show Y e for (0.5, 0.6) by the green line, for (1.35, 1.2) 
by the blue line, and (1.8, 2.4) by the maroon line. For the (1.35, 1.2) case, we had noted before 
in Figure 2, that Pp e = P£ e = 1 over the entire energy range. For the case (1.8, 2.4) there is no 
conversion as well. Hence for these case there is no flavor conversion and Y e stays constant for 
all r, given solely in terms of the 4> r 9e /4> Ve and ((E — A np ) 2 ) Pe / ((E + A„ p ) 2 )„ e ratios. For the case 
(0.5, 0.6) (cf. Figure 1) we have single split in both the neutrino and antineutrino channels. For 
this case therefore we see a variation in Y e as a function of the radius. Since we have noted that 
for cases of {4> r v , <f>p) for which there is flavor conversion due to collective effects, Y e fluctuates 
non-trivially with the radius for r ^ 400 km, therefore in what follows, we will show all results for 
r ^ 400 km. That implies that we consider only the neutrino-driven wind region henceforth. 

Figure 8 also shows the most significant aspect of flavor conversion due to collective effects. 
At r = 30 km, one is in the synchronization phase and the collective bipolar oscillations are yet to 
set in. This is the limiting case of no conversion already discussed, while by r = 400 km, they are 
complete. The maroon line (1.5, 2.0) for IH and green line (0.5, 0.6) for NH show that collective 
effects can change the value of the electron fraction from Y e < 0.5 (at r = 30 km) to Y e > 0.5 (at 
r = 400 km). Hence, we can explicitly see that these combination of values of (jf Ve and <p r De (for the 
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Figure 9: The exclusion plot consistent with Y e < 0.45 for the spectrum Gl, G3 and LL for both NH and IH. 
The allowed area is to the left of the curves. The dotted orange lines denote the Lower Limit (LL) and Upper 
Limit (UL) of </>£ and <jf D for Gl and G3, arising from the two fold uncertainty defined in Eq. (25). Similarly 
double dotted dashed brown lines denote the Lower Limit (LL) and Upper Limit (UL) for the Lawerence Livermore 
spectrum model (LL). 



respective hierarchies) will not allow r-process once collective effects are taken into account and 
hence will be ruled out if one imposes the criteria that Y e must be less than 0.5. Therefore, it is 
expected that the exclusion plot for any specific limit on Y e will change once collective effects are 
taken into account. 

In Figure 9 we show the exclusion plot in the </>y e -</>£ e plane for Y e < 0.45. The lines themselves 
correspond to the case when Y e = 0.45, while the allowed area (which gives Y e < 0.45 or \p eP / K e n > 
1.22) is to the left of the curves. We reiterate that the survival probabilities for these plots 
have been calculated at r = 400 km. We vary <p Ve and <fr Pe in the ranges (0.5,5.0) and (0.5,3.5) 
respectively. We plot exclusion curves for both IH and NH, and for all the three spectra models. 
We also plot the exclusion curves when collective oscillations are absent (WOC) for LL (thick red 
dashed) and G1/G3 (thick sky blue long dashed). These lines correspond to the no conversion 
case discussed earlier. We can see that they are almost straight lines in the 4>l e -(f) 9e , with 4> r Pe /(fi Ue = 
1.22* ((E + A np ) 2 ) Ue / ((E — A np ) 2 ) Pe , for the respective spectral model. Since the average energy as 
well as the spectral shape (cf. Eq. (17)) for both Gl and G3 are the same, the exclusion lines for 
no oscillation case for them is identical. For LL, since the ratio of ((E + A np ) 2 ) Ve / ((E — A np ) 2 ) Pe is 
smaller (0.62) than for G1/G3 (0.88), the exclusion line for no oscillation case for LL corresponds 
to smaller <fr r Pe / '(j) r v . We see from the figure that for NH the effect of collective oscillations are 
mainly in the low (f? Ue and <j) Pe region. This conforms to the observation in section 2, where we had 
shown that the split energy increased with increasing <f) r Ue and/or <p P< _. Since the flux begins to fall 
with increasing energy, the impact of collective oscillations fall for higher (jf and <j) Pe . Hence for 
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Figure 10: Exclusion plots for NH and IH for model G3. The exclusion condition is varied from Y e < 0.35 to 
0.5. The dotted sky blue lines denote the l W (Lower Limit) and l uV (Upper Limit) from the assumed two fold 
uncertainty of 0£ e and 0p c , for the spectrum model G3. The allowed region is on the left side of the exclusion 
curves. 

all the three models, the NH exclusion curves at higher relative luminosities agree with the WOC 
ones, as there is no observed collective effect there. Whereas at lower relative luminosities the 
curves deviate from the WOC ones due to the observed single splits, making the allowed region 
smaller. Again, as seen in Figure 8, the effect of collective oscillations is to increase Y e for any 
given (f) r 9e i '4> r v . Since Y e decreases with 4>l e /(f)l e , the exclusion plots shift to larger 4> r Pe /4> r Ue values 
once collective oscillations are switched on. This results in the curves shifting left in the <j) 9 -<fi^ 
plane. 

Figure 9 shows that for IH, the effect of collective oscillations can be very significant in con- 
straining the relative luminosities. This can be understood from Figures 7 ,8 and Eq 31. In Figure 
8 we see that the effect of collective oscillation is to shift Y e to a higher value. Eq. 31 shows that 
this can be balanced by increasing 0£ e compared to 4> r Ue . This results in shifting the contour plots 
to the left of the <f) 9 — <p v plane in Fig.9. Since Figures 7 and 9 show the same plane, we 

can see that for IH once collective oscillations are switched on, the only region in this plane which 
remains allowed is the one where we have double splits in the neutrino sector and no splits in the 
antineutrino sector. This is the (II, 0) zone where J z < and D z < 0. 

One understands that stronger constraint on Y e reduces the allowed parameter space. Figure 
10 shows the exclusion plots for IH and NH for the model G3 for various constraints Y e < 0.35, 
0.40, 0.45, 0.50. The left panel is for IH and the right one for NH. Since higher values of <fi r v gives 
higher Y e , as one reduces the required value of Y e , 0£ e gets more constrained. The constraint on 
0p e is relatively weak. But very low values of <ft Pe are not allowed as we have seen that lower values 
of <j) Pe increases the electron fraction. 

Some final comments are in order. We would like to reiterate here that our analysis and the 
corresponding exclusion plots are meant as a proof of principle and merely indicate the ranges 
of the allowed fluxes for which one gets neutron-rich regions for r-process in the neutrino driven 
wind. They show that all r-process calculations should take collective effects into account and 



20 



such detailed simulations can be used to extract more rigorous bounds on the initial fluxes. 

Now we discuss the evolution as a function of time and its implication for our results. Early 
times have larger luminosities and may deviate more from energy equipartition. In realistic sim- 
ulations in the cooling phase the luminosity for each species decreases with time as mentioned 
earlier but it may be reasonable to assume that the relative luminosities change very slowly with 
respect to time. With time the shock moves out, the neutrinospheres slowly fall in, the matter 
in the hot bubble and the wind driven region cools and the constituents change, first producing 
alpha-particles and then heavier nuclei. The alpha particles are strongly bound systems and their 
excitation by vjv can be neglected. The inclusion of a-particles in the matter was considered in 
[47] and the effect of the nuclear composition involving heavier nuclei on Y e was looked at in [16]. 
A similar study including spectral splits in a self-consistent manner for the time-evolved system 
need to be undertaken separately in future. 

One also realizes that there are many uncertainties that exist presently in the occurrence of the 
r-process in the supernovae. Firstly large entropy needed for the development of the r-process, as 
mentioned earlier, need to be observed consistently in all one dimensional simulations as well as in 
simulations going beyond one dimension. The hot bubble region does have the problem of having 
lower entropy [52]. This is compounded by the inability of simulations to give rise to outgoing 
shocks with the right explosion energy. Often the physics understanding comes from the use of 
'semi-analytic models' [52] that critically depend on the three quantities Y e , the entropy and the 
dynamic timescale. 

4 Summary and Conclusions 

Collective flavor oscillations driven by neutrino-neutrino interaction at the very high density region 
of core collapse supernovae control the emitted flux of neutrinos of different flavors. In the process 
one or more swaps of flavors for both neutrinos and antineutrinos take place depending on the 
initial neutrino flux and distributions. We study the phenomena of spectral splits and consequent 
flavor swaps for different models of neutrino spectrum, varying the relative luminosities of neutrinos 
and antineutrinos for both normal and inverted mass hierarchy. The effect of spectral splits is 
found to be more pronounced for inverted hierarchy and depending on the initial luminosities one 
can get single or dual splits in neutrinos and/or antineutrinos. For a specific choice of relative 
luminosity we also find a case for inverted hierarchy where the splits are not resolvable numerically 
and is akin to no spectral splits for all practical purposes. Single split patterns are also obtained 
for normal hierarchy for some choices of the luminosities. Next we consider the impact of the 
collective oscillations and the spectral splits on the electron fraction Y e , which determines if the 
environment is neutron-rich and compatible with r-process nucleosynthesis or not. The minimal 
requirement for r-process is the electron-to-nucleon ratio Y e < 0.5, but a more favorable condition 
may be Y e < 0.45/0.40. We consider the flavor evolution reduced to an effective two flavor model 
with oscillation between v e and v x and their antiparticles. The oscillation parameters are chosen 
as Am 2 = 3 x 10 -3 eV 2 and a small effective mixing angle 9 = 10~ 5 in agreement with realistic 
1 — 3 mixing. For these values of parameters the ordinary MSW resonances take place beyond 
the r-process region and hence there will be no effect of neutrino conversions on Y e due to these. 
However, the inclusion of collective effects can affect the value of Y e even for these values of mass 
and mixing parameters. 
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The electron fraction (Y e ) as a function of the radius of the core is calculated and it shows an 
oscillatory behavior in the bipolar region due to collective effects, before saturating to a constant 
value which depends on the initial luminosities and the pattern of flavor swap. Different models 
of neutrino energy distributions are used. For each of the distributions initial fluxes of different 
flavors are varied and constraints on the initial neutrino fluxes consistent with successful r-process 
nucleosynthesis are shown in exclusion plots for these initial neutrino fluxes. While a detailed 
simulation of the r-process nucleosynthesis inside the supernova might bring some changes to the 
exclusion plots, this work illustrates the fact that such exclusion plots are possible to achieve. 

The variation in the number of spectral splits with the variation in the luminosity give rise 
to different possibilities of neutrino and antineutrino spectrum at the detector. The constraints 
on luminosities obtained by ensuring r-process nucleosynthesis can provide additional inputs in 
narrowing down the possible patterns. 

Note added: While this paper was under peer review, a few papers [56, 57] have appeared which 
study the possible effects of a three flavor treatment of the collective oscillations of supernovae 
neutrinos. The differences of this with a two flavor treatment and their impact on the r-process 
nucleosynthesis need to be investigated. However from the above studies it is expected that the 
extra effect due to three flavors may be observable in IH but the changes in NH will be minor. 
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